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Abstract . Density- height profiles in the stratosphere and mesosphere 
were obtained from impact-pressure, velocity, and altitude measurements made 
on six X-15 research airplane flights. A form of the Rayleigh pitot formula 
was used for density computations. Because of pressure-instrumentation 
limitations and pressure lag, the maximum altitude for reasonably accurate 
density determination was considered to be about 65 km* Good agreement was 
obtained between temperatures calculated from faired density-height profiles 
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of two X-15 flights and temperatures measured by rocket sondes launched near 


the times of flight from the Pacific Missile Range, Point Mugu, California 
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OBTAINED WITH THE X-15 AIRPLANE 
Earl J. Montoya and Terry J. Larson 
NASA Flight Research Center, Edwards, California 

INTRODUCTION 

1 The X-15 research program provides an opportunity to obtain atmospheric 

2 data in the stratosphere and the mesosphere by virtue of both the high altitudes 

3 obtained by the X-15 airplane (Figure l) and the research instrumentation 

4 installed in support of the basic research program. Analysis showed that, with 

5 the particular instrumentation available, density could be determined more 

6 accurately than temperature or pressure. The calculation technique selected 

7 is based on^a form of the Rayleigh pitot equation similar to that used by 

8 Ainsworth , et al . [ 1961 ] for the continuum- flow regions of an Aerobee flight 

9 launched from Fort Churchill, Canada, in 1958. This method requires only 

10 impact pressure, velocity, and altitude measurements at high supersonic speeds 

11 and is, thus, particularly applicable to X-15 data, inasmuch as the X-15 

12 impact-pressure measurements are not subject to flow angularity effects. 

13 A number of density profiles were determined. Six of these profiles, from 

14 flights made over Southern California and Nevada between March 30, 19^1, and 

15 August 22, 1963, are presented and discussed in this paper. The profiles are 

16 significant in that they (l) add to the atmospheric data in the altitude range 

17 not covered by radiosonde balloons and satellite measurements and only recently 

18 being extensively explored by the Meteorological Rocket Network, and they 

19 (2) demonstrate the applicability of impact-pressure measurements from flight 

20 vehicles in the determination of atmospheric density. 
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METHODS AND MEASUREMENTS 

1 For continuum- flow conditions at supersonic speeds, atmospheric density 

2 may be obtained from a modified form of the Rayleigh pitot equation (see 

3 Appendix for derivation) 

p 3 _ p2(p _ o.460p) + 0,1 7 ^ ■■ p + 0.0525^7- » 0 

V V° 

4 where p = density, p = ambient pressure, = impact pressure, and 

5 V = velocity. The constants in this equation were evaluated, using a specific- 

6 heat ratio of 1.4 for air. Although the terms which include static pressure 

7 are relatively insignificant at the high Mach numbers (greater than 2.5) 

8 associated with the data of this paper, the preceding five-term equation was 

<4 

9 programed to obtain maximum accuracy for density determination over the entire 

10 supersonic Mach number range of the X-15. 

11 Impact-pressure measurements were obtained from the stagnation port of the 

12 X-15 airflow sensor, commonly referred to as the ball nose (Figure 2). The 

13 sensor has a pressure-nulling type of servomechanism which constantly orients 

14 the stagnation port toward the'* local airstream. Thus, the measurement of 

15 impact pressure is not subject to flow angularity effects as are pressure 

16 measurements from some so, unding rockets. 

17 Ambient-pressure measurements were obtained up to altitudes of 30 km from 

18 rawinsonde balloons released at Edwards, California, by the Air Weather Service. 

20 At the higher altitudes (for which rawinsonde data were unavailable), values 

21 from the U.S. Standard Atmosphere, 1962, were used. As previously mentioned, 

22 however, accurate determination of static pressure is not critical for the 

23 data of this paper. 

2k Velocity and altitude measurements 'were obtained from three AFMTC-Mod II 

25 radars (modified SCR-584 radars). Usually, simultaneous data were available 
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from two radars, although computations of velocity and altitude were based on 
a single- station solution. Fairings were made by subjective analysis. True 
air velocity was computed by vector addition of wind velocity and radar- 
determined velocity. The wind velocity was acquired either from AN/GMD-1 
rawinsonde measurements or from sounding-rocket measurements from the Pacific 
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Missile Range v . 

0 0 

Measurements were obtained on two of the three X-15 aircraft, the X-15-2 

and the X-15-3. On the two' density-data flights of the X-15-2, impact pressures 

* 

were not recorded directly but were obtained from a static-pressure measurement 
and a differential-pressure measurement. Static pressure was measured from the 
fuselage manifolded orifices and differential pressure from the static source 
and the stagnation port of the ball nose. These two measurements , when summed, 
gave the required impact pressure. 

Impact pressures were obtained similarly for X-15-3 data flights at alti- 
tudes below 40 km. Above this altitude, however, impact pressures were recorded 
directly by standard NACA airspeed- altitude recorders, which axe mechanical- 
optical instruments using both absolute and differential cells. Table 1 lists 
the altitude ranges of the cells. 


. LIMITATIONS AND ACCURACY 

For the X-15 flight conditions, real-gas effects on the validity of the 
Rayleigh pitot equation are insignificant. Therefore, no attempts were made 
to correct for deviations from the ideal-gas relations or for minor ionization 
effects that occur in the speed range where the measurements were made. The 
effects of the variation of specific-heat ratio for air (caloric effects) were 
also considered. By using a specific heat ratio of 1.4, it was found that 
density errors were less than 1 per cent for the .X-15 flight conditions 
encountered; thus, no caloric corrections were made. 
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According to Kane and Maslach [1950] and Sherman [ 19531 > corrections to 
impact pressures are necessary when the Rayleigh pitot equation is used for 

3 compressible flow conditions associated with viscous forces that are signifi- 

4 cant in comparison to the inertia or pressure forces is the fluid. Data from 

5 these studies indicate that corrections for viscous effects are appropriate 

6 when free- stream Reynolds numbers are less than about 300 for impact-pressure 

7 installations similar to those of the X-15 airflow sensor. For the 1.27-cm- 

8 diameter orifice of the sensor and the X-15 flight velocities encountered, 

9 Reynolds numbers below 300 are limited to altitudes greater than about 60 km. 

10 Since the basic measurement errors are, evidently, considerably larger than the 

11 indicated viscous corrections for the pertinent flight conditions, no corrections 

12 for viscosity have been applied. Because of the probable inaccuracies in the 

13 pressure data and in space positioning, the density profiles presented are 

14 smoothed density- variation curves and do not portray small-scale details of the 

15 atmospheric soundings . 

16 Of primary importance in obtaining the derived density is the measurement 

17 of impact pressure. Meaningful density data have been determined to altitudes 

18 of approximately 65 km by means of pressure recorders routinely installed on 

19 research aircraft at the NASA Flight Research Center. The cells of these 

20 recorders are accurate from 0.25 per cent to about 0.50 per cent of full-scale 

21 range. For the extreme flight conditions encountered with the X-15, the latter 

22 value was used for the estimation of density errors shown -in Table 2. The 

23 recorded pressures for the portions of flights that had simultaneous recordings 

24 of either two absolute or two differential cells of different ranges indicate 

25 disagreements significantly less than 0.50 per cent of the higher range cells. 

26 Because of space limitations behind the impact port of the ball nose, 

27 longer tubing (60 cm) than desired was required between the impact-pressure 

28 port and the recorder. As a result of the large impact opening (1.27 cm), the 
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large-bore tubing (O .75 cm minimum), and the small instrument volumes, however, 
the associated pressure lag was relatively small up to an altitude of 6 l km. 

Inasmuch as the individual errors listed in Table 2 are maximum errors 
and are independent, a 3 ° error is included to represent the overall maximum 
probable density error. Although it is not likely that the individual errors 
are additive, the fact that the maximum error due to impact-pressure measurement 
is about an order of magnitude larger than any of the other errors lessens the 
probability of significant -compensation at the higher levels of measurement. 

In addition to the density errors, it is estimated that maximum altitude errors 
are 0-3 km. 
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PRESENTATION AND ANALYSIS OF DATA 
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The geographical .location of the area (shaded) over which the measurements 
presented herein were obtained is shown in Figure 3* . A time history of the 
pertinent flight quantities (impact pressure p^, velocity V, and altitude h) 
is presented in Figure 4 to indicate the approximate range and variation of 
these quantities for a typical X-15 flight. 

The variation of density with altitude, as determined by measurements 
obtained on six X-15 flights, is shown in Figures 5 through 10. The data for 

•t 

the two X-15-2 flights (Figures 5 and 6) are not considered to be as accurate 
at altitudes above about 40 km as those of the X-15-3 flights (Figures 7 to 10) 
because of the less sensitive pressure cells on the X-15-2. A distinction is 
made between the ascent and descent portions of the flights to aid in inter- 
preting the data. Also shown are density values obtained from rawinsonde 
measurements- (at the lower altitudes) and the U.S. Standard Atmosphere, 1962. 

At altitudes below 6 l km for which lag effects should be insignificant, 
the difference between ascent and descent points. measures the repeatability of 
the data for a given altitude and flight. Below 65 km, the average scatter is 
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28 
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approximately ±5 per cent. Prom Table 2 it is evident that scatter between 
ascent and descent points for altitudes less than approximately 50 km reflects, 
primarily, velocity and altitude errors rather than pressure errors. At these 
altitudes, therefore, the density-height fairings have been made to correspond 
more closely to the points calculated from the radar data considered to be the 
most accurate. Also, more weight was given to the rawinsonde densities (where 
available), inasmuch as they are regarded as generally the most accurate of the 
data presented. 

At altitudes in excess of 65 km, the ascent and descent points diverge 
with altitude. This divergence is caused primarily by lag in the pressure- 
measuring system, as indicated by the consistently higher ascent than descent 
impact pressures at comparable altitudes. Because of this lag effect, as well 
as the known limitations of the impact-pressure measurements and viscous effects 
at high altitudes, meaningful density data for any of the flights cannot be 
expected at altitudes above 65 km. 

The validity of the faired density- height profiles at altitudes where 
meaningful data can be expected was indirectly confirmed by data obtained from 
two Areas rocketsondes launched from the Pacific Missile Range near the times 
of two X-15-3 flights. Temperatures calculated from the density- height profiles 
of Figures 9 and 10 are compared in Figures 11 and 12, respectively, with 
temperatures measured by the Areas 10-mil bead thermistor. 

Densities were converted to temperatures by using the- equation 



where M is -the molecular weight, R is the gas constant, 
ation of gravity, and dh is the differential of altitude, 
refer to the end points of the considered altitude interval. 


g is the acceler- 
The subscripts 
By using the 
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temperature given by the U.S. Standard Atmosphere, 1962, at 6l km and the 
densities from the curves of Figures 9 and 10 for this altitude, temperatures 
were calculated for progressively lower altitudes in intervals of 0.6l km for 
each flight. A constant value of molecular weight (M = 28.9 for air) was used. 
As pointed out by Ainsworth et al . Cl96l3 , an error in the initial choice of 
temperature becomes relatively insignificant after integration over a number 
of kilometers of altitude. 

The Areas temperature data were corrected for radiation effects as sug- 
gested by Wagner [ 1963 ]* The small differences between Areas and X-15 tempera- 
tures (Figures 11 and 12) indicate good agreement, considering the possible 
errors in the Areas temperature measurements (±2 per cent with Wagner's correc- 
tions) and the ^horizontal separation (maximum of about 375 km) between the two 
sets of data. 

In Figures 5 through 10, the density-height profiles consistently indicate 
higher densities than the standard values above altitudes from about 35 km to 
45 km. Figure 10 indicates a possible crossover point at altitudes in excess 
of approximately^60^km7^ As shown, the differences between the flight data and 
the Standard Atmosphere generally vary consistently from about 5 to 7 per cent 
over large altitude ranges. This consistency adds further validity to the 
density-height fairings and the usefulness of the measurements made with the . 
X-15 airplane. 


APPENDIX 


DERIVATION OF ATMOSPHERIC DENSITY FROM RAYLEIGH PITOT EQUATION 

The relationship between stagnation pressure p^. measured behind a normal 
shock, free-stream static pressure p, and free-stream supersonic Mach number M 
is given by the Rayleigh pitot equation 


lii 


f*(7 + 1) 

7-1 — =7 

■y 7-1 M 2 

1 ^ 7 ' 1 ^1 

L 2 J 

7 M 

27 M 2 


1 

7-1 


where 7 is the specific-heat ratio for the gas. Proceeding as outlined by 
Newell [l953]) — may be substituted for a 2 (speed of sound) so that 

p 

-P pV 2 * 

Vr = . Substituting this expression into the Rayleigh pitot equation and 

expanding 


H = BpV 2 
■p 7P 


1 . _e_ + -seL + 7 . (, 2 ?. 

2PV 2 8pV^ 48p^V°^ 


where 


B = 


(7 + !)j 


Zli 

7-1 


._z_ 

7-1 


which may be expressed as the cubic 


p3-^( Pt .^) + 2^ + i^llp3 
bv 2 ' ' 6 v km 


0 


By choosing 7 = lA for air, the derived equation becomes 


p 3 - i#p2(p t - 0.4605P ) ♦ + „ 0 

TT V V 
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TABLE 1 . Pressure-Measurement Cell Ranges 


Range 

Airplane 

Altitude, 

km 

Absolute cell, 
mm Hg 

Differential cell, 
mm Hg 

X-15-2 

16— 26 

0—790 

0—940 


26—38 

0—30 

0—940 


38-66 

0—30 

0—90 

X-15-3 

18—4.0 

0—800 

0 

& 

1 

0 


40—74 

0—40 

Not required 
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TABLE 2. Estimation of Density Errors 





Density errors, % 



Altitude, 

km 

Velocity 

P 

Pt 

Lag 

7 

Max. probable 
errort 

30-5 

5-2 

Negligible 

2.4 

Negligible 

0.75 

5-8 

45-6 

7-5 

2.0 

2. 5(6. 8)* 

Negligible 

0.75 

8.2(10.4)* 

61.0 

2.6 

2-5. 

11.2(28.0)* 

1.0 

0.75 

11.8(28.3)* 


*X-15-2 only. 

t3a (root sum square) . 
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Fig. 1. X-15 airplane in the landing configuration. 






Fig. 3* Geographical area over which X-15 density measurements were made. 









Fig. 7 . Atmospheric-density-height profile determined from X-15-3 flight 

on May 2, 1963- 








Fig. 10. Atmospheric-density- height profile determined from X-15-3 

flight on August 22 , 1963* 
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Data 

O X -15 flight 
D Corrected Areas 
O Edwards rawinsonde 

Date of flight: August 22, 1963 

X-15 launch time: 1606 Z 

Areas launch time: 1622 Z 

Rawinsonde release time: l600 Z 



T, °K * 

Fig. 12. Comparison of temperatures calculated from density-height profile 
of figure 10 with Areas rocketsonde temperatures. 
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WAVELENGTH, microns 



(a) Typical oscillogram from time-of-f light scanning spectrometer. 

(b) Calibration curve of time-of-flight scanning spectrometer. 

(c) Reduced experimental spectra. 


Fig. 4.- Steps in determination experimental spectra. 
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Fig. 6.- Typical shadowgraph of hypersonic projectile and assumed shock shape. 
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